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We consider a nonlinear version of the Landau-Zener problem, focusing on photoassociation of 
a Bose-Einstein condensate as a specific example. Contrary to the exponential rate dependence 
obtained for the linear problem, a series expansion technique indicates that, when the resonance is 
crossed slowly, the probability for failure of adiabaticity is directly proportional to the rate at which 
the resonance is crossed. 
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When external conditions force the diabatic energies of 
two interacting quantum states cross, under conditions 
of complete adiabaticity the system makes a transition 
from one bare energy eigenstate to the other. However, 
"complete adiabaticity" applies when the levels cross "in- 
finitely slowly". The realistic case when only a finite 
time is spent near resonance makes the Landau- Zcncr 
(LZ) problem |||§, one of the handful of key paradigms 
in quantum mechanics. On the other hand, the Gross- 
Pitaevskii equation and its variants, nonlinear versions 
of standard quantum mechanics, have proven highly suc- 
cessful in the theoretical description of alkali- vapor Bosc- 
Einstein condensate [0]. Nonlinear quantum mechanics 
gives up the superposition principle and hence the linear 
vector space structure. How much of our intuitive un- 
derstanding of quantum mechanics survives becomes an 
issue. 

The nonlinear version of the LZ problem that is the fo- 
cus of this Letter is something we have come across in our 
studies of photoassociation of a Bose-Einstein con- 
densate (BEC). However, the same question will come 
up in an attempt to control the scattering length Jq] of 
an atomic BEC by means of a Feshbach resonance 0, in 
second-harmonic generation in nonlinear optics |||, and 
generally in field theories where the nonlinearity is of the 
form ift cfxf). For concreteness we couch our discussion in 
terms of photoassociation, and thus consider processes 
in which two atoms and a laser photon join to make a 
molecule. It is known from a second-quantized numerical 
solution that by sweeping the frequency of the photoas- 
sociating laser, in the adiabatic limit it is possible to con- 
vert an atomic BEC entirely into a molecular BEC [Q. 
The question is, what are the nonlinear consequences if 
the frequency is not swept "infinitely slowly"? While 
we actually know of no publication addressing this ques- 



tion [p 10 1 , we have found that the same applies also to 
the semiclassical version of the problem, when the atomic 
and molecular condensates are described by a nonlinear 
two-state system. 

The answer to the above question is outlined as follows. 
First the essentials of the linear LZ problem are reviewed, 
which serves both as an introduction to the formalism 
and a foundation for our series expansion. Turning to 
the nonlinear LZ model, we use photoassociation as an 
example to review the dressed state picture and illustrate 
the physics of adiabatic following. In particular, a pic- 
ture in terms of crossing of stationary states still applies 
and predicts near-adiabaticity, but the details are qual- 
itatively different from the case of linear quantum me- 
chanics. The first-order Heisenberg equations of motion 
for the semi-classical atom-molecule amplitudes are then 
converted into a third-order differential equation for the 
molecular probability, whereby a perturbative expansion 
allows for what amounts to a Green function solution. 
All told, we find that, while in the linear case the prob- 
ability for no transition is an exponential function of the 
speed at which the resonance is crossed, at slow sweep 
rates of the laser frequency the probability that an atom 
does not make a transition to the molecular condensate 
is directly proportional to the rate. Before closing, we es- 
timate explicit numbers for two systems, 23 Na and 87 Rb, 
that could serve to test our prediction. 

The linear LZ problem is specified by the equations 



ia = Q(3, i$ = ila + 8(3 . 



(1) 



Here a and [3 are the probability amplitudes for the 
two bare or diabatic states, Vl is the (real) coupling be- 
tween them, and the detuning 5 — £t stands for the 
time dependent frequency difference between the two 
states. For any fixed 8 the system (fy) has two nonde- 



1 



generate dressed |TT[ or adiabatic states [a,0\ with the 
property that the time dependent solution is of the form 
[a{t),3{t)\ — e~ let [a,Po] for some quasienergy he. In 
particular, one of these states is such that [a, 3) — ► [1,0] 
when S — ► -co, and evolves continuously to [0, 1] as 
8 — > oo. 

Suppose that the system with the time dependent de- 
tuning 8 = £i starts out in the state with a = 1 when 
i — > — oo, then it also start out in this dressed state. As 
the detuning is swept "infinitely slowly" to +oo, by the 
adiabatic theorem the system emerges at time t = oo 
in the same dressed state, or with the probability am- 
plitudes that satisfy \a\ — 0, \3\ = 1. More accurately, 
when the detuning is not swept infinitely slowly, £ > 0, 
adiabaticity tends to break down especially in the neigh- 
borhood of 8 = where the two quasienergies have their 
closest approach. The time evolution of the probabil- 
ity for state 3, P(t) = \0(t)\ 2 . often written Q using 
parabolic cylinder functions |l2j| , can more conveniently 
be expressed in terms of the Kummer and Tricomi con- 
fluent hypergeometric functions \Fi and U ]l2| : 

P L z{\t) = \-e-* x \U(-i\± -^t 2 )| 2 , t<0, (2a) 



P L z{\t) = l-e^ A | 



+ t/(-a i,-i^ 2 )| , t>o, 



(2b) 



where T is the gamma function. At t = +oo we have the 
Landau-Zener result 



where 



P? z (\)=P LZ (\oo) = l-e 



x- Q2 



-4ttA 



(3) 



(4) 



is the conventional LZ parameter. 

Let us next turn to the nonlinear LZ problem, for 
which we develop a series expansion solution. Modern 
mathematical physics is of course replete with series ex- 
pansion techniques. For example, scries expansions have 
been used to solve an ion in a Paul trap the 3- 

state Potts |Q and Heisenberg anti-ferromagnetic 
models, scattering jl6| and gauge theories j[7), molecular 
adsorption JlJ], black hole physics as well as quan- 
tum anharmonic [^0|,^l) and double well potentials 
Moreover, a series expansion approach has been applied 
to systems with cubic [||| and quartic [^| nonlinearites, 
which are directly related to the nonlinear version of the 
Landau-Zener problem considered herein. 

The equations to solve are a semiclassical approxi- 
mation to the photoassociation problem p3], which de- 
scribes atomic and molecular condensates not as boson 
fields but as classical fields. We have 



ia = —= a* 3, 
y/2 



ip = — — a 
y/2 



8(3. 



(5) 



Physically, a and 3 are the probability amplitudes that 
an atom is part of the atomic or the molecular conden- 
sate. Equations (|^), in fact, preserve the normalization 
|a| 2 + |/3| 2 = 1. The Rabi frequency for coherent pho- 
toassociation, fi, may be adjusted by varying the inten- 
sity of the driving laser field. The detuning 5, again 
swept linearly as S — measures the mismatch between 
the laser frequency and the frequency in a transition in 
which two atoms are converted into a molecule. It is 
controlled by tuning the photoassociating laser. Transi- 
tions to noncondensate modes |25|-|2"7|] are neglected un- 
der the assumption that the condensate coupling satisfies 
f2 <C hp 2 / 3 /m P,p|,p6[| , where p is the initial atom density 
and m is the atomic mass. 

For a fixed detuning the system (j|) has dressed states 
as well H , though they behave quite differently from the 
dressed states of the linear system ([!]). They are of the 
form [a(t), 8(t)] — [e~ v *<5, e~ 2v */3], where Tip is conven- 
tionally referred to as chemical potential. Depending on 
the value of the detuning, there may be as many as three 
essentially different dressed states, along with a multi- 
tude of additional states that can be generated using the 
symmetry transformations of Eqs. (JsJ) . A trivial state M 
exists for all detunings. It has all atoms as molecules, 
(3 = 1, and the chemical potential equals p = ^8. A 
second state B exists for S/il < \/2 and has a = — 1 at 
8 = — oo. For this state /3 — > 1 and p — > ^8 as 8/il — > 2, 
so that at <5/f2 = V2 the states M and B continuously 
merge. The third state F exists for 8 /ft > —\/2 and 
turns into all atoms, a — * 1, in the limit 8 — > oo . 

Suppose now that the system starts as all atoms at 
a large negative detuning, i.e., in state B, and that the 
detuning is swept slowly through the resonance. One 
expects adiabatic following, so that the system should 
stay in state B and eventually turn into state M, into 
molecules. This actually is the result both in the cor- 
responding quantum case addressed in Refs. pL[|, and 
from our numerical trials with Eqs. (^|). However, the 
point of resonance where the degree of the transfer of 
population is mostly determined occurs where the chem- 
ical potentials are equal, at S/Q = \/2, which is also 
precisely where the dressed states merge. The situation 
is mathematically much different from the linear LZ case. 

It is therefore gratifying that we have been able to 
carry out an analysis analogous to the LZ case in our 
nonlinear system. The (quite involved) technical details 
will be reported elsewhere, here we only present the key 
idea. Thus, it is possible to combine from Eqs. (|5j> a 
third-order differential equation for the probability for 
molecules, P = |/3| 2 , which reads 
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Here the primes refer to derivatives with respect to the 
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transformed time variable z = £t' 
equation for the linear problem ( 
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V(\P) = 0, (7) 



where V refers to the rule for forming the left-hand side 
of the differential equations The solutions to the latter 
are, of course, known; Plz(z) of Eqs. (|J) is the one with 
the appropriate initial conditions at t = — oo. 
Now, Eq. (ra) may be rearranged to read 



V{X,2P) = 3A 



^4P'P 



p2 

^2 



(8) 



In the limit with A < 1 the LZ transition probability is 
small, and so is presumably the solution to Eq. @. In 
Eq. (^) the nonlinear terms, the same ones that make 
the right-hand side of Eq. (g) , are compared with terms 
that are much larger. The nonlinearity makes a "small" 
perturbation. We are thus lead to formulate the Ansatz 



P{z) = \P L z(\z) + eP {1) {z)+e 2 P {2 Xz) + 



(9) 



with e = Pf^ z may be thought of as the small parameter 
in a perturbative expansion. 

The calculations make use of what is essentially the 
Green's function for the linear differential operator acting 
on P in Eq. (0). It turns out that a formal expansion of 
the type (j^) may be found, and for a small enough A it 
even converges to the desired solution of Eq. (^). The 
leading terms in the expansion at infinite time are 
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Equation (^) may also be rearranged as an equation 
for the atomic state probability \a\ 2 = R = 1 — P, and 
in this form reads 



£>(-|A,P) = 3A - 



2R'R P 2 -l/3 
_+ 2z 2 



(11) 



In the limit A> 1 very few atoms will remain, so one is 
tempted to attempt an Ansatz of the form 

R(z) = [1 - Plz(±A,z)]+ e rW(z) + e 2 R^{z) + ... . 

(12) 

This time the formally small parameter is e — R(oo) 
1 — Pj^ z . Perhaps surprisingly, this Ansatz works, too. 
By employing essentially the same mathematics as in the 
case of the expansion (^), one finds a formal expansion 
that may be shown to converge for large enough A. The 
leading terms give 
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(13) 



The most notable qualitative difference from the linear 
case occurs in the limit of slow sweep of the frequency, 
A ^> 1. The probability for no transition, |a(oo)| 2 
I — P(oo), behaves in the linear case like exp(— iirX), but 
in the nonlinear case as (37rA) _1 . 

We have solved Eqs. (||) also numerically using several 
methods. An example is shown in Fig. |l](a), where we 
plot the atomic probability R = 1 — P = \a\ 2 as a func- 
tion of the running detuning S = t;t. In this example we 
initially fix a negative detuning that lies outside the range 
of the graph, and start the system in the corresponding 
dressed state B. We then slowly increase the rate of the 
sweep of the detuning until it reaches £ = 0. 1 Q 2 , still 
far to the left of the detuning axis in Fig. [l](a). Upon 
arrival to the range plotted in Fig. 0(a), the detuning is 
swept at this constant rate £, and the occupation prob- 
ability is within 10~ 8 of the occupation probability for 
the dressed state B. At this point adiabaticity still pre- 
vails. Especially around 5 /SI ~ a/2, though, adiabatic- 
ity breaks down. To demonstrate, we show in Fig.|l(b) 
the difference between the actual occupation probability 
of the atomic state and the prediction from the dressed 
state B. Typical of both linear and nonlinear numerical 
solutions, oscillations arise because the system does not 
perfectly follow the ground state, and filtering out these 
oscillations accurately is a major tour de force. Nonethe- 
less, it is easy to see the A -1 dependence of the residual 
atomic probability on the sweep rate numerically. 

While these aspects of the computations are somewhat 
trivial, computations also provide new insights. The lin- 
ear LZ problem is completely symmetric in the two states. 
Whether the system starts in state a or /3, the probabil- 
ity for the transition to the other state is the same. This 
cannot hold in the nonlinear case; if a — at some time, 
then by Eqs. ([bJ) it will remain so forevermore. Also, 
given that the states B and M merge at 5 = \/2 fi, one 
might expect that there is some preponderance to tran- 
sitions to a molecular state. This does not seem to be 
the case. In Fig. |l|(c) we resort to the same scheme as 
in Fig. |l|(a), except that we start the system in the su- 
perposition with a = (3 = . The probabilities oscillate 
violently around S — 0, but eventually settle in the neigh- 
borhood of half and half again. 

We now briefly consider explicit numbers for 23 Na 
and 87 Rb condensates. First is the validity of the 
two- mode model (m, which again hinges on the con- 



dition ft <g; hp 2 / 3 }m. The condensate coupling is ||: 
fl = [(A 3 p)(/// )] 1/2 , where 2ttA is the wavelength of the 
photoassociating light and Iq is the chracteristic inten- 
sity of a photoassociation transition with binding energy 
~ 1 cm -1 . Given the 23 Na ( 87 Rb) characteristic intensity 
I Q = 0.47 W/cm 2 (I = 0.07 W/cm 2 ) @, and a typical 
BEC density p — 5 x 10 14 cm -3 , the two-mode model is 
valid for / < 50mW/cm 2 (/ < 10mW/cm 2 ). The lin- 
ear rate dependence for the transition probability should 
manifest for slow sweep rates, i.e., £ ~ 0.01 SI 2 . For in- 
tensities satisfying the above relation, and typical BEC 
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densities, this translates into a 23 Na ( 87 Rb) detuning- 
sweep rate £ ~ 200 x 2ttHz (£ - 50 x 2ttHz). 

In sum, we have studied a nonlinear variant of the 
Landau-Zener problem. Though our specific example 
was about photoassociation of a Bose-Einstein conden- 
sate, the problem is generic in classical and bosonic field 
theories with a cubic nonlinearity. While the basic adia- 
baticity argument still works, the structure of the adia- 
batic or "dressed" states in the nonlinear system is much 
different from its counterpart in standard linear quantum 
mechanics. When the levels cross slowly, the probability 
for the failure of adiabaticity turns out to scale linearly 
with the rate at which the levels cross, as opposed to the 
exponential behavior of the linear problem. This predic- 
tion could be tested with existing 23 Na and 87 Rb exper- 
imental setups. 
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FIG. 1. (a) Probability of remaining in the atomic con- 
densate R — 1 — P as a function of running detuning 5 = ££, 
with £ = 0.01 Q 2 . The system is started at a large negative 
time in the dressed eigenstate B, and arrives adiabatically 
from the left to the range of 5 shown in the figure, (b) Dif- 
ference between the atomic probability R and the adiabatic 
probability Rb of the dressed state B for the same data as 
in panel (a). We set Rb = for 8 > \/2fl, when the state 
B has merged with the all-molecules state M. (c) Same as 
panel (a), except that the system was started in a superposi- 
tion with a = (3 = 1/ y/2. 
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